Quantum Transport with Spin Dephasing: 
A Nonequlibrium Green's Function Approach 
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A quantum transport model incorporating spin scattering processes is presented using the non- 
equilibrium Green's function (NEGF) formalism within the self-consistent Born approximation. 
This model offers a unified approach by capturing the spin-flip scattering and the quantum effects 
simultaneously. A numerical implementation of the model is illustrated for magnetic tunnel junction 
devices with embedded magnetic impurity layers. The results are compared with experimental data, 
revealing the underlying physics of the coherent and incoherent transport regimes. It is shown that 
small variations in magnetic impurity spin- states/concentrations could cause large deviations in 
junction magnetoresistances. 

PACS numbers: 72.10.-d, 72.25.-b, 72.25.Rb, 71.70.Gm, 73.43.Qt 
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I. INTRODUCTION 



Quantum transport in spintronic devices is currently 
a topic of great interest. Most of the theoretical work 
reported so far has been based on the Landauer ap- 
proach [l[ assuming coherent transport, although a few 
authors have included incoherent processes through aver- 
aging over a large ensemble of disordered configurations 
0) However, it is not straightforward to include 

dissipative interactions in such approaches. The non- 
equlibrium Green's function (NEGF) formalism provides 
a natural framework for describing quantum transport 
in the presence of incoherent and dissipative processes. 
Here, a numerical implementation of the NEGF formal- 
ism with spin-flip scattering is presented. For magnetic 
tunnel junctions (MTJs) with embedded magnetic impu- 
rity layers, this model is able to capture and explain three 
distinctive experimental features reported in the litera- 
ture 0,01 regarding the dependence of the junction 
magnetoresistances (JMRs) on (1) barrier thickness, (2) 
barrier height and (3) the number of magnetic impurities. 
The model is quite general and can be used to analyse 
and design a variety of spintronic devices beyond the 1-D 
geometry explored in this article. 



This article is organized as follows. In the view of ped- 
agogical clarity, a heurisic presentation of the NEGF for- 
malism with spin dephasing mechanisms is given Sec [TT1 
followed by a numerical implementation of the model in 
Sec IIIII Initially (Scc lIII A"]) the definitions of device char- 
acteristics are presented for impurity free MTJs together 
with device parameters benchmarked against experimen- 
tal measurements. How to incorporate the spin exchange 
scattering mechanisms into the electron transport calcu- 
lations is shown (Sec IIIIBj) . and the model is applied to 
MT J devices with magnetic impurity layers (Sec MI Cjl . 
Theoretical estimates and experimental measurements 
are compared as well in this section (Sec IIII CI) , while 
a summary of the results is given in Sec IIVI 



II. MODEL DESCRIPTION 

NEGF Method: The problem is partitioned into chan- 
nel and contact regions as illustrated in Fig. Q] [9] . Com- 
ponents of the partitioned device can be classified in four 
categories: 

(i) Channel properties are defined by the Hamiltonian 
matrix [H] including the applied bias potential. 

(ii) Contacts are included through self-energy matrices 
[SrJ/pij] whose anti-hermitian component: 
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FIG. 1: A schematic illustration of device partitioning in 
NEGF formalism. Magnetization direction of the drain is 
defined relative to the source(A# = #r— #l). 



describes the broadening due to the coupling to the con- 
tact. The corresponding inscattering/outscattering ma- 
trices are defined as: 

XT,r (E) =fo(E- » LtR ) T L , R (E) , (2a) 



^Tr (E) = [l-fo(E- (il,r)] T l , r (E) , (2b) 

where f (E - (j, L}R ) = 1/1 + cxp [(E - ^ L ,R)/k B T] is 
the Fermi function for the related contact. 

(iii) Electron- electron interactions are incorporated 
through the mean field electrostatic potential matrix [U] . 

(iv) Incoherent scattering processes in the chan- 
nel region are described by in/out-scattering matrices 



[S'cf] /[Sg ut ]. Broadening due to scattering is given by: 
T S (E) = [E™ (E) + Eg* (E)] , (3) 

from which the self-energy matrix is obtained through a 
Hilbert transform: 
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Eqs. (Q][3]) are the boundary conditions that drive the 
coupled NEGF equations [Eqs. (JSKH])], where Green's 
function is defined as: 

G= [EI-H -U -Vl-Vr-VsT 1 , (5) 

with the spectral function (analogous to density states): 

A = i[G-G^] = G n + G P , (6) 

where [G™]/[G P ] refer the electron/hole correlation func- 
tions (whose diagonal elements are the electron/hole den- 
sity): 



G n, P = q 



i in, out 
J R 



Gl (7) 



The in/out-scattering matrices [E™]/[Eg u *] are re- 
lated to the electron/hole correlation functions [G n ]/[G P ] 
through: 



G"* (r,r';E T h^)d(hu;) 



- / E ^ r '-> M] „ s/ (r, r'; 25) d (M 



(8) 



Here the spin indices (crfe, 07) refer to the (2x2) block di- 
agonal elements of the on-site electron/hole correlation 
function which is related through the [Z) n ]/[Z) P ] tensors 
to the (<7i,<Xj-) spin components of the (2x2) block di- 
agonal of the in/out-scattering function. This term can 
be interpreted as in the following. The first part de- 
scribes the process of spin-flip transitions (subscript sf) 
due to the spin-exchange scatterings in the channel re- 
gion. The second part denotes the contributions of the 
spin-conserving exchange scatterings (subscript nsf for 
"no spin-flip") in the channel region. Both of the con- 
tributing parts are previously shown by Appelbaum[To| 
using a similar treatment. Here the [£>™]/[Z? p ] are fourth- 
order scattering tensors, describing the spatial corre- 
lation and the energy spectrum of the underlying mi- 



croscopic spin-dephasing scattering mechanisms. These 
scattering tensors can be obtained from the spin scatter- 
ing hamiltonian: 



ff i^ = E fl . J^-R^Z-Si, (9) 

where r/Rj are the spatial coordinates and a / Sj are the 
spin operators for the channel electron/ (j-th) magnetic- 
impurity. For point like exchange scattering processes, 
the scattering tensor for the first term in Eq. © cor- 
responding to the spin-flip transitions is given by (see 
Appendix) : 



3 



[D">p (r, r>\ fkj)] at = S(r - r') £ ^N^) 





ITT) 


III) ITI) 


lit) 


(TTI 





Fu.dS (w =F w g ) 


0" 


(III 


F d ,u8 (UJ ± LOq) 








(Til 











(ITI 








0. 



(10) 



Also, the no-spin flip component in Eq. ((SJ is given by: 



(ct^jI I 

[D">p (r, r'; ^)] nsf = 6(r r>) £ | J 5 ^, (w - w,)] 



(TTI 
(III 
(Til 
(ITI 



ITT) 


III) 


ITI) 


IIT) 


l 








" 





1 














-1 














-1 



(11) 



Current is calculated from the self-consistent solution 
of the above equations for any terminal "i": 



oo 




trace ([£*" (E) A (E)] - [T< (E) G n (E)])dE). 



— OO 

(12) 

The general solution scheme without going into the de- 
tails can be summarized as follows . The matrices listed 
under the categories (i) and (ii) are fixed at the outset of 
any calculations. While the [U], [£™' ou *] and [Eg] ma- 
trices under the charging and scattering categories (iii) 
and (iv) depend on the correlation and spectral func- 
tions requiring an iterative self- consisted solution of the 
NEGF Equations [Eqs. (3j|5])]. One important thing to 
note is that for the numerical implementation presented 
in Sec |ni] we do not compute the charging potential [U] 
self-consistently with the charge. The change in tunnel 
barriers is neglected and assumed not to influence the 
electrostatic potential. This allows one to focus on the 
dephasing due to the spin-flip interactions. 



the impurities have an internal degree of freedom, such 
as fluctuating internal spin states in the case of mag- 
netic impurities. Electron scatterings from such impuri- 
ties randomize the electron spin states which can not be 
simply incorporated through the device Hamiltonian. In- 
stead scattering self energy matrices are needed for this 
type of phase-relaxing scattering proccesses. An imple- 
mentation of this self-energy matrix treatment will be 
discussed in this part of the paper for electron-impurity 
exchange scattering processes in MTJs. But first, we'll 
discuss MTJ fundamentals and device characteristics in 
the absence of magnetic impurities (the coherent regime). 



A. MTJs: Coherent Regime 

MTJs devices considered here consist of a tunnel- 
ing barrier (A10 x ) sandwiched between two ferromag- 
nets (Co) with different magnetic coercivities enabling 
independent manipulation of contact magnetization di- 
rections (Fig. [2]). Single band tight-binding approxi- 



III. APPLICATION: MTJs WITH MAGNETIC 
IMPURITIES 

In the presence of "rigid" scatterers such as impuri- 
ties and defects, electron transport is considered coher- 
ent as the phase relationships between different paths 
are time independent. Accordingly, any static devia- 
tions from perfect crystallinity leading to the scattering 
of the electrons from one state to another can be incor- 
porated into the transport problem through the device 
Hamiltonian H. However the situtation is different when 




FIG. 2: Energy band diagram for the model MTJs considered 
here. 
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mation is adopted [llj with an effective electron mass 
(to* = TO e ) in the tunneling region and the ferromag- 
netic contacts. Accordingly for constant effective mass 
throughout the device, transverse modes can be included 
using 2-D integrated Fermi functions /2D [E z — ^l,r) — 
(m*k B T/2nh 2 ) In [1 + cxp {ji LtR - E z /k B T)] in Eq. © 
instead of numerically summing parallel k-componcnts. 
The Green's function of the device is simply defined as: 

G(E Z ) = {E Z I- H- E L -E B ) -1 , (13) 



without any self-consistent solution requirements where 
H is Hamiltonian of the isolated system, and E^/E^ 
are the self-energies due to the source/drain contacts. 
Here the matrices have twice the size of the channel re- 
gion in the corresponding presentation due to the elec- 
tron spin states. Accordingly, in real space representa- 
tion for a discrete lattice whose points are located at 
x = 3 a i j being an integer (j = 1---N), the matrix 
(E Z I — H — El — E/j) can be expressed as: 



|1> 



\N-l) 



\N) 



E Z I-H-X, 



E fl 



<1| 
(2| 

(N-l\ 

(N\ 



EJ-cci-Y.L (3 

(3+ E z I-a 2 







E Z I -a N -i (3 

13+ E z I-a N 



(14) 



where a n is a 2x2 on-site matrix: 
2t + U n 



El 

c,n 









- 2t + U n 



(15) 



and j3 = —tl is a 2x2 site-coupling matrix with t = 
h 2 /2ma 2 and I = ( J ^ ). The left contact self-energy 
matrix is nonzero only for the first 2x2 block: 



. te ik\a 





_ te ik l L a 



(16) 



where E z = E ] c ' 1 + U L + 2t (l - cos k^ l aj . For the right 
contact only the last block is non-zero: 



E fl (N,N,E Z 



-te lk R a 






-te ik R a 



f , (17) 



where E z = E^ 1 +U R + 2t( 1 — cosfcj^a) with §t being 
the unitary rotation operator defined as: 



ft(A0) = 



cos(A6»/2) sin(A6»/2) 
-sin(A6»/2) cos(A6»/2) 



(18) 



for two contacts with magnetization directions differing 
by an angle Ad. 

A theoretical analysis of MTJ devices in the absence of 
magnetic impurity layers is presented and compared with 
the experimental data 0, 0] for varying tunneling barrier 
heights and thicknesses. The parameters used here for 
the generic ferromagnetic contacts are the Fermi energy 



E F = 2.2 eV and the exchange field A = 1.45 eV [llj . 
The tunneling region potential barrier [J7b arr ] is parame- 
terized within the band gaps quoted from the literature 
0, G3 1 while the charging potential [U] is neglected due 
to the pure tunneling nature of the transport. 

Coherent tunneling regime features are obtained by 
benchmarking the experimental measurements made in 
impurity free tunneling oxide MTJs at small bias volt- 
ages. Referring to If/Iaf as the current values for the 
parallel/antiparallel magnetizations (A6 = 0/A6 = it) of 
the ferromagnetic contacts, the JMR is defined as: 



JMR = (I f -Iaf)/If. 



(19) 



The dependence of the JMRs on the thickness and 
the height of the tunneling barriers is shown in Fig. [3Ja) 
with an energy resolved analysis [Fig. [3[b)] for different 
barrier thicknesses (0.7 — 1.4 — 2.1 nm). JMR values are 
shown to be improving with increasing barrier heights for 



all barrier thicknesses, a theoretically predicted | 7|, |8 l , 14 1 
and experimentally observed 0, [11 Gil, E HHH ISf, 
[22l . I23I [24[ feature in MTJs. The barrier heights obtained 
here may differ from those reported in literature [l5l [7j| 
[g|ll El M, HI, S3, [H S3 based on empirical models 

Experiments and theoretical calculations observe de- 
terioration of JMRs with increasing barrier thicknesses 
[Fig. GJa)]. Whereas an energy resolved theoretical anal- 
ysis shows that energy by energy junction magnetoresis- 
tances defined as: 



JMR{E Z ) = (I F (E Z ) - I AF {E Z ))/I F {E Z ), 



(20) 



remain unchanged [Fig. (2Kb)]. This initially counter in- 
tuitive observation can be understood by considering the 
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FIG. 3: For impurity free MTJs, (a)thickness dependence of 
JMRs for different barrier heights are shown in comparison 
with experimental measurements [l5l . [l6l Il8l . fTsL |20| . [2ll . 
E3, H3] while an (b)energy resolved analysis of JMR(E Z ) 
(left-axis) and normalized w(E z ) (right-axis) distributions is 
also presented for a device with a tunneling barrier height of 
1.6 eV. For MTJs with impurity layers, (c) variation of JMRs 
for varying barrier thicknesses and interactions strengths are 
shown together with (d)an energy resolved analysis. Normal- 
ized JMRs are proven to be thickness independent as dis- 
played in the inset. 



redistribution of tunneling electron densities over ener- 
gies with changing tunneling barrier thicknesses. Defin- 
ing w(E z ) as a measure of the contributing weight of 
the JMR(E Z ), one can show that experimentally mea- 
sured J MR is a weighted integral of JMR(E z )s over E z 
energies: 



J MR = J oj(E z ) J MR (E z ) dE z , 



(21) 



where lu (E z ) = Ip (E z )/Ip is the energy resolved spin- 
continuum current component (weighting function). In 
Eq. (f2Tj) , independently from the barrier thicknesses 
JMR(E Z ) ratios are constant (solid line in Fig. [31(b)), 
while the normalized w(E z ) distributions shifts to- 
wards higher energies with increasing barrier thicknesses 
(dashed lines in Fig.[3Kb)). Hence JAIRs, an integral of 
the multiplication of the w (E z ) distributions with the en- 
ergy resolved J MR (E z )s, decreases with increasing bar- 
rier thicknesses (Eq. (jT9j)) . 



B. Adding Spin Exchange Scattering 

Spin exchange scattering processes are responsible 
from the incoherent nature of the tunneling transport 
for the model devices considered here. Through this elas- 
tic scattering process, the state of nothing else seem to 
change if the electron and the impurity spins are con- 
sidered as a composite system. Nevertheless, it is an 
incoherent process since the state of the impurity has 
changed. What makes this process incoherent are the 
external forces forcing the impurity spins into local equi- 
librium. The incoherent nature of the scattering lies in 
the "information erasure" of the surrounding through 
the forces forcing the impurity spins into an unpolarizcd 
(%50 up, %50 down) spin distribution. These external 
forces in a closely packed impurity layer can be magnetic 
dipole-dipole interactions among the magnetic impuri- 
ties or spin relaxation processes coupled with phononic 
excitations. Nevertheless the physical origin of the equi- 
librium restoring processes is not of our interest (at least 
from tunneling electron's point of view) assuming equilib- 
rium restoring processes are fast enough to maintain the 
impurity spins in a thermal equilibrium state. Accord- 
ingly, NEGF formalism incorporates spin dephasing ef- 
fects of the environment into the electron transport prob- 
lem through a boundary condition (the spin-scattering 
self-energy). As discussed in Sec HH coupling between 
the number of available electrons/holes ([G™]/[G P ]) at 
a state and the in/out-flow ( [Ef] /[££"*]) to/from that 
state is related through the fourth order scattering tensor 
[D n ]/[DP] in Eq. ©. 

For the model systems considered here magnetic im- 
purity spin states are degenerate (ftw q — 0) allowing only 
elastic spin-flip transitions. Spin-conserving scattering 
processes are also neglected due to their minor effect on 
the JMRs. Accordingly, the spin scattering tensor rela- 
tionship given in Eqs. (|8I10I and fTTj) will simplify to: 



(r, r; E) 



sf 



Glf^r-E). (22) 



[Z?"]/[Z? P ] scattering tensors relate the electron/hole 
[G™]/[G P ] correlation matrices with the (2x2) block diag- 
onal elements of [Eg 1 ] /[Eg 1 *] in/out-scattering matrices 
of form: 
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through the tensor relationship (shown below in matrix 
format) for the corresponding lattice site "j" with mag- 
netic impurities: 



^SiTT J, 



J S;TT 

\T8iU ). 



-^in.out N 



(^in,out\ 



j j 



j j 



J Z N, 





Fd, u 





F u ,d 






(°tf) 



where iVV 



(24) 

impurities and 



is the number of magnetic 
Fu/Fd represents fractions of spin- up/spin-down impu- 
rities for an uncorrelated ensemble (F u + Fd = 1). 

This tensor relationship can be understood heuristi- 
cally from elementary arguments. The in/out-scattering 
into spin-up component is proportional to the density of 
the spin-down electrons/holes times the number of spin- 
up impurities, N[F U : 



J S\TT 



= J'NtR 



(25) 



while preserving the normalized u{E z ) carrier distribu- 
tions [Fig.EUd)]. 

Further analysis shows that "normalized" JM Rs 
deteriorates with increasing spin-dephasing strenghts 
(J 2 Ni) independently from the tunneling barrier heights 
[Fig. Ufa)]. This general trend can be shown by mapping 
the "normalized" JMRs into a single universal curve us- 
ing a tunneling barrier height dependent scaling constant 
c(U bar r) (inset in Fig. Ufa)). 

This allows us to choose a particular value of barrier 
height [Ubarr — 1.6eV] and adjust a single parameter J 
to fit our NEGF calculations [Fig. BJb-d)] with experi- 
mental measurements obtained from (5-doped MTJs 
Submonolayer impurity thicknesses given in the measure- 
ments are converted into number of impurities using ma- 
terial concentrations of Pd/Ni/Co impurities and device 
cross sections (6 x 10~ 4 cm 2 ) [5j. 

Close fitting to the experimental data are observed at 
77K [Fig. @Jb-d)] using physically reasonable exchange 
coupling constants of J = 2.63/iey/2.68/zeVyi.41/ieV for 
devices with Pd/Ni/Co impurities [2(|. However, experi- 
mentally observed temperature dependence of normalized 



Similarly, the in/out-scattering into spin-down com- 
ponent is proportional to the density of the spin-up elec- 
trons/holes times the number of spin-down impurities, 
NiF d : 



JJ 



jj 



(26) 



C. Incoherent Regime: Results 

The incoherent tunneling regime device characteristics 
in the presence of magnetic impurities is studied for a 
fixed barrier height Ubarr = 1-6 eV [Fig. [3] (c)] with 
changing barrier thicknesses and electron-impurity spin 
exchange interactions (J 2 N] = — 6 eV). Nonlinear de- 
creasing JMRs with increasing spin-exchange interac- 
tions are observed at all barrier thicknesses due to the 
mixing of independent spin-channels |5, 6] while the nor- 
malized JMRs are proven to be thickness independent 
(inset). This observation is attributed to the elastic na- 
ture of the spin exchange interactions yielding a total 
drop in JMR(E Z ) values at all E z energies in Eq. 




0.2 0.4 0.6 
10" 12 N. 



FIG. 4: (a) Normalized JMRs deteriorates with increasing 
spin-dephasing strengths independently from the tunneling 
barrier heights. This general trend can be scaled to a single 
universal curve (inset). Experimental data taken at UK and 
300A' is compared with theoretical analysis in the presence of 
(b) Pd, (c) Ni and (d) Co magnetic impurities with increasing 
impurity concentrations. 



7 



JM R ratios can not be accounted for by our model calcu- 
lations. Broadenings of the electrode Fermi distributions 
due to changing temperatures from 77K to 300K seem 
to yield variations in normalized JMR ratios within a 
linewidth. As a result, different J exchange couplings 
are used in order to match the experimental data taken 
at 300K. 

For Pd and Ni doped MTJs, relatively small variations 
in J exchange couplings are needed (J300/J77 — 1-16 for 
Pd and J300/J77 = 1-19 f° r Ni) in order to match the 
experimental data at 300K. These small temperature de- 
pendences could be due to the presence of some secondary 
mechanisms not included in our calculations. One such 
mechanism reported in literature includes the presence of 
impurity-assisted conductance contribution through the 
defects (possibly created by the inclusion of magnetic im- 
purities within the barrier) which is known to be strongly 
temperature dependent [27j. In fact, the contribution of 
the impurity-assistant conductance is proportional with 
the impurity concentrations in accordance with the ex- 
perimental measurements. 

Another interesting feature observed in calculations 
[Fig. |l{b-c)] is the comparable J exchange couplings for 
the Pd and Ni impurities at temperatures 77K and 300K. 
Accordingly, a possible estimate of the impurity spin 
states may be made by considering the most commonly 
encountered oxidation states of the Pd and Ni impuri- 
ties. Closed-shell elemental Pd is only known to be in a 
magnetic oxidized state of S=l in octahedral oxygen co- 
ordination according to the Hund's rules [28| . Similarly, 
we may attribute the comparable J couplings for Ni im- 
purities due to the S=l spin state of the Ni +2 which is 
known to be a frequently observed ionized state in oxygen 
environment [29| . 

On the contrary, for Co doped MTJs, there's a clear 
distinction for normalized JMR ratios at different tem- 
peratures [Fig. 0Jd)] which can not be justified by the 
presence of secondary mechanisms. Fitting these large 
deviations require large variations in J exchange coupling 
parameters [^300/^77 = 1.73], We propose this to be a 
result of thermally driven low- spin/high- spin phase tran- 
sition [33, [3l[ . This is a credible argument since the oxi- 
dation state of the cobalt atoms can be in Co +2 (S = 3/2, 
high-spin) or Co (S = 0, low-spin) state or partially in 
both of the states depending on the oxidation environ- 
ment. Such thermally driven low-spin/high-spin phase 
transitions for metal-oxides have been predicted by theo- 
retical calculations and observed in experimental studies 
[3fl l3~fl | . These phase transitions have not been discussed 
in MTJs community in connection with possible scatter- 
ing factors determining the temperature dependence of 
JMRs. Although from the available experimental data 
it's not possible to make a decisive conclusion in this 
direction, given the nonlinear dependence of JMRs on 
magnetic impurity states in our calculations, we believe 
that it's important to point out this possibility here. 



IV. SUMMARY 

Summary: A NEGF based quantum transport model 
incorporating spin-flip scattering processes within the 
self-consistent Born approximation is presented. Spin- 
flip scattering and quantum effects are simultaneously 
captured. Spin scattering operators are derived for the 
specific case of electron-impurity spin-exchange interac- 
tions and the formalism is applied to spin dependent elec- 
tron transport in MTJs with magnetic impurity layers. 
The theory is benchmarked against experimental data in- 
volving both coherent and incoherent transport regimes. 
JM Rs are shown to decrease both with barrier thickness 
and with spin-flip scattering but our unified treatment 
clearly brings out the difference in the underlying physics 
[Fig. [3]. Our numerical results show that both barrier 
height and the exchange interaction constant J can be 
subsumed into a single parameter that can explain a va- 
riety of experiments [Fig. EJa)] . Nonlinear dependence 
of JMR ratios on the impurity concentrations are shown 
[Fig. BJb-d)] in difference with earlier linear estimates 
from empirical models Accordingly, small differences 
in spin-states/concentrations of magnetic impurities are 
shown to cause large deviations in JMRs. Interesting 
similarities and differences among devices having Pd, Ni 
and Co impurities are pointed out, which could be sig- 
natures of the spin states of oxidized Pd and Ni impuri- 
ties and low-spin/high-spin phase transitions for oxidized 
cobalt impurities. 



Acknowledgments 

This work was supported by the MARCO focus cen- 
ter for Materials, Structure and Devices and the NSF 
Network for Computational Nanotechnology. Part of 
this work was carried out at the Jet Propulsion Labo- 
ratory under a contract with the National Aeronautics 
and Space Administration NASA funded by ARDA. 



APPENDIX: SPIN DEPHASING SELF ENERGY 

In the following, the scattering tensors stated in 
Eqs. (jlOlllip of the main paper are obtained starting from 
standard expressions obtained in the self-consistent Born 
approximation from the NEGF formalism. Here we start 
from the formulation in Ref. d (see Sections 10.4 and 
A. 4) which represents a generalization of the earlier treat- 
ments [HI, H3, HU, Hfl . For spatially localized scatterers 
we have: 



Hi maj (t) = (a l \H I (t)\a j 



(A.2b) 



D n aiaj . t(Tk(Jl (t,t') = ( Y T'twe (*')) 
= ( Y (t) \s p ) (sp\H\ m<T . (0 | Sq )\ (A.la) 



4^(0 = ^l^ T (0k> (A-2c) 



rt /v 



= ( E (^i^^ w i s «) < s °i (0 M ) ( A - lb ) 



where \s a ) and |s^) are impurity spin subspace states and 
Hj is the interaction Hamiltonian [Eq. (J9j> ] defined within 
the channel electron spin subspace as: 



Hl- ai a h (t) = <<7i|J?j(t)|<Tk> 



<„ ) .(t') = <^|ff](t')ki> (A.2d) 
For an uncorrelated impurity spin ensemble with: 

P = Y Ws « I s ") ( s «l = E Ws « ^1 ( A - 3 ) 



averaging in Eqs. (|A.la - lA~Tb ) can be done through a 
weighted summation of spin scattering rates of magnetic 
(A. 2a) impurities: 



D^..^ (t,f) = Y, w ^ (* a \H!, aiah (t) \s p ) (s \ H\ maj (f) \s a ) 

= Y MHiw* (*) \p]H\ maj (f) \s a ) = tr ( P H\. ai(T . (f ) F J;CTi(Tfc (t) 



(A.4a) 



Dl ia ^ kai (t,t') = Y W *° ( S /3l#Wi (*) l«a> (^l^ka, (*') M 

s a ,sp 

= Y IfiAHiw, (*) M#U«r» (0 M = *r (p^^ (*') (*) 



(A.4b) 



The trace tr (pA) for any operator A is independent of 
representation. Accordingly, [D n ]/[D P ] scattering ten- 
sors can be evaluated using any convenient basis for mag- 
netic impurity spin states. 

Through Jordan- Wigner transformation, single spins 
can be thought as an empty or singly occupied fermion 
state: 



IT)^a+|0), 



l!> = |o>, 



(A.5a) 



(A.5b) 



with creation/annihilation operators for the channel elec- 
trons: 



a+(t) =a + (t) -- 
a(t) = a~{t) = 



e iUet 



-iu e t o 



(A.6a) 
(A.6b) 



For degenerate electron spin states (fuo e = 0), there's no 
time dependence as such <v (t) — > <r , a (t) — > a. 

Pauli spin matrices are related to creation/annihilation 
operators through a x = (a' + a)/2, o~ y = (at + a)/2i and 
a~ z = a^a — 1/2. Accordingly, interaction Hamiltonian 

Hi (f,t) = J 8 (r — R) a ■ S (t) can be expressed as: 



9 



Hi 



(r, t) = JS (f- i?) ^aS+ (t) + i a tS_ (t) + (eta - Sz (t) 



(A.7) 



r 



where S is the spin operator for the localized magnetic Substituting the interaction Hamiltonian from 

impurity. Eq. (TA~7j) into Eqs. fOk - TOb ) will yield: 



K^) -> 

(o»C.j| 4 



(o^,^ (*) = 



(TTI 
(III 
(til 
(III 



(TTI 
(III 
(Til 
(ITI 



ITT) 

S z (f ) S z (t) 
5- (f ) 5+ (t) 
S_ (f ) S z (t) 
Sz (f ) 5+ (t) 

ITT) 

Sz (*') Sz (t) 
5+ (f ) S_ (t) 
Sz (f ) S_ (t) 
5+ (f) S z (i) 



III) 

5+ (f ) S_ (t) 
Sz (f ) S z (i) 
-Sz («0 S_ (i) 
-S + (*') S z (t) 

III) 

S- (t')S + (t) 
Sz (?) Sz (t) 
-S_ (f ) Sz (t) 

-s z (Os+(t) 
I 





ITI) 




lit) 


s+ 


(*') Sz 


(t) 


Sz (to s_ (t) 1 


Sz 


(t')S + 


(t) 


-S_ (f ) Sz (t) 


Sz 


(f) Sz 


(t) 


s_ (tf)si (t) 


s+ 


(t')S + 


(t) 


-Sz (tO Sz (t) . 




III) 




Ill) 


Sz 


(f) s + 


(t) 


S_ (f ) Sz (t) 1 


s+ 


(tO Sz 


(t) 


-Sz (to S_ (t) 


Sz 


(tO Sz 


(t) 


s- (tos- (t) 


s+ 


(tO 5+ 


(t) 


-Sz (tO Sz (t) . 



(A.8a) 



(A.8b) 



The localized magnetic impurity spin-operators can be 
written in its diagonalized impurity spin-subspace as: 



S 4 







(A.9a) 



with iVj being total number of impurities, the desired 
quantities [D n ]/[D P ] can be obtained by evaluating the 
expectation values of the operators in Eqs. (|A.8a[) and 
()A.8b[) . Here the only non-zero elements are: 



S_ = d = 



S z = dU 







1 

-1 



(A.9b) 



(A.9c) 



with uj g = AEj/H where AEj is the energy difference 
between spin-up and spin-down states for the localized 
magnetic impurities. 

For a given impurity density matrix of the form 
(Fu+Fd=l): 



tr( P Sz(t')Sz(t)) = -, (A. 11a) 

tr( P S+(t')S-(t))=F u e- iu «( t - t '), (A.llb) 



tr ( P S_ (t')S+(t))=F d e l ^( t - t '). 



(A.llc) 



P = Nj(u 9 ) 



F u 
F d 



Finally, for a given impurity density matrix 
(A. 10) [Eq (jA~T0)) ]. [D n ]/[D p ] tensors are obtained as: 



\<Jkcri) — > 

(o-i<Tj\ i 



D n (t,t') = J2J 2 N I (u 



(TTI 
(III 
(Til 
(ITI 



F d e 



ITT) 


III) 


ITI) 


lit) 


1/4 












1/4 














-1/4 














-1/4 J 



(A.12a) 
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1 \ 




ITT) 


1 1 1 \ 

III) 


III) 


1 1 +\ 

III) 
















(TTI 


1/4 








DP (t,t') = E^/K) 


(III 
(Til 


Fu g-iw,(t-t') 




1/4 





-1/4 








(III 











-1/4 J 



(A.12b) 



It's convenient to work with the Fourier transformed 
functions as such (t — t') — » hco: 



e -i Uq (t-t>) e -n\t-t'\/H ^ 6 r HuJ _ faj 



(A.13) 



where rj being a positive infinitesimal. With 
Fourier transforming Eqs. (|A.12allA.12b[) will simplify to 



Eqs. (jlOHlip [35l ] . For the calculations reported in this 
article, diagonal elements not leading to spin-dephasing 
are omitted due to their negligible effect on J MR ratios. 
In this case [£)™]/[Z) p ] scattering tensors simplifies to a 
form (Eq. (|24p) which can be understood from simple 
common-sense arguments (Eqs. (|25|) and ((26)) ). 
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